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Abstract 



i We develop a field theoretical approach to the cold interstellar medium 

i (ISM) and large structure of the universe. We show that a non-relativistic 

self-gravitating gas in thermal equilibrium with variable number of atoms 
or fragments is exactly equivalent to a field theory of a single scalar field 
. (/^{x) with exponential self-interaction. We analyze this field theory per- 

Q>^ ! turbatively and non-perturbatively through the renormalization group ap- 

proach. We show scaling behaviour (critical) for a continuous range of the 
Qh' temperature and of the other physical parameters. We derive in this frame- 

Q ■ work the scaling relation M{R) ~ R'^" for the mass on a region of size R, 

and Av ~ R'^ for the velocity dispersion where q = ^{dn — 1). For the 
, density-density correlations we find a power-law behaviour for large distances 

~ \fi — r2p'^^^^- The fractal dimension dn turns to be related with the 
■ critical exponent of the correlation length by dfj = The renormaliza- 

^ • tion group approach for a single component scalar field in three dimensions 

^ states that the long-distance critical behaviour may be governed by the (non- 

perturbative) Ising fixed point. The Ising values of the scaling exponents 
are = 0.631..., dn = 1.585... and q = 0.293.... Mean field theory yields 
for the scaling exponents v = 1/2, dn = 2 and q = 1/2. Both the Ising 
and the mean field values are compatible with the present ISM observational 
data: 1.4 < dn < 2, 0.3 < q < 0.6 . As typical in critical phenomena, the 
scaling behaviour and critical exponents of the ISM can be obtained without 
dwelling into the dynamical (time-dependent) behaviour. We develop a field 
theoretical approach to the galaxy distribution. We consider a gas of self- 
gravitating masses on the FRW background, in quasi-thermal equilibrium. 
We show that it exhibits scaling behaviour by renormalization group meth- 
ods. The galaxy correlations are first computed assuming homogeneity for 
very large scales and then without assuming homogeneity. In the first case we 
find ^(r) =< p{fo)p{fQ + r) > / < p —1 ~ r~^ , with 7 = 2. In the second 
case we find D{r) =< p(ro)p(ro -|-r) > ~ r~^ with F = 1. While the universe 
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becomes more and more homogeneous at large scales, statistical analysis of 
galaxy catalogs have revealed a fractal structure at small-scales (A < WOh~^ 
Mpc), with a fractal dimension D = 1.5 — 2 (Sylos Labini et al 1996). 

We study the thermodynamics of a self-gravitating system with the theory 
of critical phenomena and finite-size scaling and show that gravity provides 
a dynamical mechanism to produce this fractal structure. Only a limited, 
(although large), range of scales is involved, between a short-distance cut-off 
below which other physics intervene, and a large-distance cut-off, where the 
thermodynamic equilibrium is not satisfied. The galaxy ensemble can be con- 
sidered at critical conditions, with large density fluctuations developping at 
any scale. From the theory of critical phenomena, we derive the two indepen- 
dent critical exponents i' and r] and predict the fractal dimension D = 1/v 
to be either 1.585 or 2, depending on whether the long-range behaviour is 
governed by the Ising or the mean field fixed points, respectively. Both set 
of values are compatible with present observations. In addition, we predict 
the scaling behaviour of the gravitational potential to be r~2(^+'?). That is, 
^-0.5 £qj, mean field or r~^'^^^ for the Ising fixed point. The theory allows to 
compute the three and higher density correlators without any assumption or 
Ansatz. We find that the connected A^-points density scales as r^^^ when 
ri >> Tj, 2 < i < N . There are no free parameters in this theory. 



I. THE SELF-INTERACTING GRAVITATIONAL GAS 

We review recent work [|l|-|^ on the statistical properties of a self-interacting gravita- 
tional gas in thermal equilibrium. We discuss two relevant astrophysical applications of 
the self-interacting gravitational gas: the cold interstellar medium (ISM) and the galaxy 
distributions. 

This review is organized as follows. In section II we summarize the main properties 
of the ISM and our results for it, in section III we develop the field theory approach to 
the gravitational gas. A short distance cutoff is naturally present here and prevents zero 
distance gravitational collapse singularities (which would be unphysical in the present case). 
Here, the cutoff theory is physically meaningful. The gravitational gas is also treated in a 
D-dimensional space. In section IV we study the scaling behaviour and thermal fluctuations 
both in perturbation theory and non-perturbatively (renormalization group approach), = 
fJ-Tgff acts as the dimensionless coupling constant for the non-linear fluctuations of the field 
0. We show that the non-linear fluctuations of the field make the theory to scale ( critical 
behaviour) for a continuous range of values of the parameters of the theory. Thus, changing 
the parameters keeps the theory at criticality. The renormalization group analysis made 
in section IV confirm such results. External forces to the ISM gas like stars are shown not 
to affect the scaling behaviour of the gas. That is, the scaling exponents q, du are solely 
governed by fixed points and hence, they are stable under perturbations. 

In section V we discuss the galaxy distributions starting from the its observational status, 
in section VI we clarify the definition of the galaxy correlators and then in section VII we 
develop the dynamical equations in the comoving frame. We apply our field-theory approach 
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in section VIII. We derive in this framework two and three point correlators deriving the 
scahng exponents. Discussion and remarks are presented in section IX. 

II. THE INTERSTELLAR MEDIUM: INTRODUCTION AND RESULTS 

The interstellar medium (ISM) is a gas essentially formed by atomic (HI) and molecular 
{H2) hydrogen, distributed in cold (T ~ 5 — 50K) clouds, in a very inhomogeneous and 
fragmented structure. These clouds are confined in the galactic plane and in particular 
along the spiral arms. They are distributed in a hierarchy of structures, of observed masses 
from 1 Mq to 10^ Mq. The morphology and kinematics of these structures are traced by 
radio astronomical observations of the HI hyperfine line at the wavelength of 21cm, and 
of the rotational lines of the CO molecule (the fundamental line being at 2.6mm in wave- 
length), and many other less abundant molecules. Structures have been measured directly 
in emission from O.Olpc to lOOpc, and there is some evidence in VLBI (very long based 
interferometry) HI absorption of structures as low as 10"*^ pc = 20 AU (3 10^^ cm). The 
mean density of structures is roughly inversely proportional to their sizes, and vary between 
10 and 10^ atoms /cm? (significantly above the mean density of the ISM which is about 
0.1 atoms/ cm^ or 1.6 10^^^ g/cm^ ). Observations of the ISM revealed remarkable relations 
between the mass, the radius and velocity dispersion of the various regions, as first noticed 
by Larson [Q, and since then confirmed by many other independent observations (see for 
example ref. |^). From a compilation of well established samples of data for many different 
types of molecular clouds of maximum linear dimension (size) R, total mass M and internal 
velocity dispersion At; in each region: 

M{R) ~ R'^» , ^v^ R\ (2.1) 

over a large range of cloud sizes, with 10~^ — 10~^ pc < R < 100 pc, 

1.4 < ci// < 2, 0.3 < g < 0.6 . (2.2) 

These scaling relations indicate a hierarchical structure for the molecular clouds which is 
independent of the scale over the above cited range; above 100 pc in size, corresponding to 
giant molecular clouds, larger structures will be destroyed by galactic shear. 

These relations appear to be universal, the exponents dn, Q are almost constant over all 
scales of the Galaxy, and whatever be the observed molecule or element. These properties 
of interstellar cold gas are supported first at all from observations (and for many different 
tracers of cloud structures: dark globules using ^^CO, since the more abundant isotopic 
species ^^CO is highly optically thick, dark cloud cores using HCN or CS* as density tracers, 
giant molecular clouds using ^^CO, HI to trace more diffuse gas, and even cold dust emission 
in the far-infrared). Nearby molecular clouds are observed to be fragmented and self-similar 
in projection over a range of scales and densities of at least 10^, and perhaps up to 10^. 

The physical origin as well as the interpretation of the scaling relations (|2.1|) are not 
theoretically understood. The theoretical derivation of these relations has been the subject 
of many proposals and controversial discussions. It is not our aim here to account for all 
the proposed models of the ISM and we refer the reader to refs. for a review. 
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The physics of the ISM is complex, especially when we consider the violent perturbations 
brought by star formation. Energy is then poured into the ISM either mechanically through 
supernovae explosions, stellar winds, bipolar gas flows, etc.. or radiatively through star 
light, heating or ionising the medium, directly or through heated dust. Relative velocities 
between the various fragments of the ISM exceed their internal thermal speeds, shock fronts 
develop and are highly dissipative; radiative cooling is very efficient, so that globally the ISM 
might be considered isothermal on large-scales. Whatever the diversity of the processes, the 
universality of the scaling relations suggests a common mechanism underlying the physics. 

We propose that self-gravity is the main force at the origin of the structures, that can be 
perturbed locally by heating sources. Observations are compatible with virialised structures 
at all scales. Moreover, it has been suggested that the molecular clouds ensemble is in 
isothermal equilibrium with the cosmic background radiation at T ~ ?)K in the outer parts 
of galaxies, devoid of any star and heating sources [0]. This colder isothermal medium might 
represent the ideal frame to understand the role of self-gravity in shaping the hierarchical 
structures. Our aim is to show that the scaling laws obtained are then quite stable to 
perturbations. 

Till now, no theoretical derivation of the scaling laws eq.( p.l|) has been provided in which 
the values of the exponents are obtained from the theory (and not just taken from outside 
or as a starting input or hypothesis). 

The aim of our work fl]-^ is to develop a theory of the cold ISM. A first step in this goal 
is to provide a theoretical derivation of the scaling laws eq. (|2.1| ), in which the values of the 
exponents dn-, Q are obtained from the theory. For this purpose, we will implement for the 



ISM the powerful tool of field theory and the Wilson's approach to critical phenomena |]T5 

We consider a gas of non-relativistic atoms interacting with each other through Newto- 
nian gravity and which are in thermal equilibrium at temperature T. We work in the grand 
canonical ensemble, allowing for a variable number of particles N. 

Then, we show that this system is exactly equivalent to a field theory of a single scalar 
field (pix) with exponential interaction. We express the grand canonical partition function 
Z as 



Z 



J J I?0e-^[^(-)] , (2.3) 



where 



Sim ^Tf^ fd'x 



T^ff = An^ , ^^ = ^lzGmy^VT^ (2.4) 

m stands for the mass of the atoms and z for the fugacity. We show that in the 0-field 
language, the particle density expresses as 

< p{f) >= --^ < VV(r) >=7^ < e^^"^ > . (2.5) 

where < . . . > means functional average over 0(.) with statistical weight e"'^'''^*^ '". Density 
correlators are written as 
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C(fi,r2) = < p(fi)p(f2) > - < p(fi) >< p(r2) > 



4 



T-I 2 



(2.6) 

'-eff ^ 

The 0-field defined by eqs. (|2.3|) - (|2.4| ) has remarkable properties under scale transformations 

X ^ xx = \x , 

where A is an arbitrary real number. For any solution 0(x) of the stationary point equations, 

VV(^) + /i' e<^(^') = , (2.7) 
there is a family of dilated solutions of the same equation (p.7|), given by 

(f)x{x) =(j){Xx) +\ogX^ . 
In addition, S[(j)x{.)] = A^"^ S[(j){.)]. 



We study the field theory ( p.3|) -( p7^) both perturbatively and non-perturbatively. 
The computation of the thermal fluctuations through the evaluation of the functional 
integral eq.( |2.3| ) is quite non-trivial. We use the scaling property as a guiding principle. In 



order to built a perturbation theory in the dimensionless coupling g = y fiT^ff we look for 
stationary points of eq. (|2.4| ). We compute the density correlator eq. (p.6|) to leading order in 
g. For large distances it behaves as 

C{n, f.) 2 r - 12 + O ( \n - r.r) . (2.8) 

32 TT^ |ri ~ ^ ' 

We analyze further this theory with the renormalization group approach. Such non- 
perturbative approach is the more powerful framework to derive scaling behaviours in field 
theory p|-0]. 

We show that the mass contained in a region of size R scales as 



<M{R)>=m J < e'^^^^ > d^x c:::^ mV a + m ^ 



+ ... 



a 

and the mass fluctuation, (AM(_R))^ =< > — < M >^, scales as 

AM{R) ~ R'^" . 

Here u is the correlation length critical exponent for the 0-theory ( p^.3| ) and a and K are 
constants. Moreover, 

K 1 

< p(r) >= ma + m ; r~~^ for r of order ~ R . (2.9) 

^ ^ 47rz/(l - a) ^ ^ 

The scaling exponent p can be identified with the inverse Haussdorf (fractal) dimension 
dff of the system 

dn = - . 

V 
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In this way, M ~ R^" according to the usual definition of fractal dimensions 
In Fourier space, 







<p{k)>r^k--. (2.10) 



From the renormalization group analysis, the density-density correlators (|2.6|) result to 

be, 

C(fi,r2) ~ |ri -fsl^-^ (2.11) 

Computing the average gravitational potential energy and using the virial theorem yields 
for the velocity dispersion. 

This gives a new scahng relation between the exponents du and q 



2 

The perturbative calculation ( p.8|) yields the mean field value for v [|l^. That is. 



i^=\ , dH = 2 and q = ^. (2.12) 

We find scaling behaviour in the 0-theory for a continuum set of values of /i^ and Teff. 
The renormalization group transformation amounts to replace the parameters /i^ and Tg// 
in j3 H and S'[0(.)] by the effective ones at the scale L in question. 

The renormalization group approach applied to a single component scalar field in three 
space dimensions indicates that the long distance critical behaviour is governed by the (non- 
perturbative) Ising fixed point |T5|-|l7|. Very probably, there are no further fixed points pO 



The scaling exponents associated to the Ising fixed point are 

1^ = 0.631... , diy = 1.585... and g = 0.293... . (2.13) 



Both the mean field (|2.12| ) and the Ising (|2.13|) numerical values are compatible with the 



present observational values ( |2.1| ) - ( |2.2| ) 



The theory presented here also predicts a power-law behaviour for the two-points ISM 
density correlation functionl (see eq. (|2.11|) , 2du — 6 = —2.830 . . ., for the Ising fixed point 
and 2du — 6 = —2 for the mean field exponents), that should be compared with observa- 
tions. Previous attempts to derive correlation functions from observations were not entirely 
conclusive, because of lack of dynamical range |]2^, but much more extended maps of the 
ISM could be available soon to test our theory. In addition, we predict an independent 
exponent for the gravitational potential correlations (~ r"^""^, where r/zsmg = 0.037. . . and 



Vmean field = |IT6|), which could bc chcckcd through gravitational lenses observations in 
front of quasars. 



The mass parameter [see eq. ( p.4| )] in the (/)-theory turns to coincide at the tree level 
with the inverse of the Jeans length 
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J 



(2.14) 



We find that in the scaling domain the Jeans distance dj grows a.s < dj >^ R. This 
shows that the Jeans distance scales with the size of the system and therefore the instability 
is present for all sizes R. Had dj being of order larger than R, the Jeans instability would 
be absent. 

The gravitational gas in thermal equilibrium explains quantitatively the observed scal- 
ing laws in the ISM. This fact does not exclude turbulent phenomena in the ISM. Fluid 
flows (including turbulent regimes) are probably relevant in the dynamics (time dependent 
processes) of the ISM. As usual in critical phenomena |T^|T^, the equilibrium scaling laws 
can be understood for the ISM without dwelling with the dynamics. A further step in the 
study of the ISM will be to include the dynamical (time dependent) description within the 
field theory approach presented in this review. 

If the ISM is considered as a flow, the Reynolds number Rcjsm on scales L ~ lOOpc 
has a very high value of the order of 10^. This led to the suggestion that the ISM (and the 
universe in general) could be modelled as a turbulent flow (Larson first observed 
that the exponent in the power-law relation for the velocity dispersion is not greatly different 
from the Kolmogorov value 1/3 for subsonic turbulence). 

It must be noticed that the turbulence hypothesis for the ISM is based on the comparison 
of the ISM with the results known for incompressible flows. However, the physical conditions 
in the ISM are very different from those of incompressible flows in the laboratory. (And 
the study of ISM turbulence needs more complete and enlarged investigation than those 
performed until now based in the concepts of flow turbulence in the laboratory). Besides 
the facts that the ISM exhibits large density fluctuations on all scales, and the observed 
fluctuations are highly supersonic, (thus the ISM can not viewed as an 'incompressible' and 
'subsonic' flow), and besides other differences, an essential feature to point out is that the 
long-range self-gravity interaction present in the ISM is completely absent in the studies of 
flow turbulence. In any case, in a satisfactory theory of the ISM, it should be possible to 
extract the behaviours of the ISM (be turbulent or whatever) from the theory as a result, 
instead to be introduced as a starting input or hypothesis. 

III. FIELD THEORY APPROACH TO THE GRAVITATIONAL GAS 

Let us consider a gas of non-relativistic atoms with mass m interacting only through 
Newtonian gravity and which are in thermal equilibrium at temperature T = f3~^. We shall 
work in the grand canonical ensemble, allowing for a variable number of particles A^. 

The grand partition function of the system can be written as 




(3.1) 



where 
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N 2 1 

i<i7^<Nh-(ij\ 



(3.2) 



G is Newton's constant and z is the fugacity. 

The integrals over the momenta pi, (1 < / < A^) can be performed exphcitly in eg. (pJ] 
using 



(2^ 



m \ 
2^ 



3/2 



We thus find, 



OO 1 



N=0 



m 



3/2' 



N 



•' 1=1 



(3.3) 



We proceed now to recast this many-body problem into a field theoretical form 
n|,[12|Jl|,|22|. 



Let us define the density 



TV 



(3.4) 



such that, we can rewrite the potential energy in eq. (|3.3|) as 
-pGm — — = -pG'm 

2 l<ljLj<N W ~ 2 J\x-y\>a 



l\x-y\>a \X - y\ 



p{x)p{y) . 



(3.5) 



The cutoff a in the r.h.s. is introduced in order to avoid self-interacting divergent terms. 
However, such divergent terms would contribute to Z by an infinite multiplicative factor 
that can be factored out. 
By using 

—An 5{x — y) , 



\x - y\ 



and partial integration we can now represent the exponent of the potential energy eq.( p.5|) 

(3.6) 



functional integral [O 

j3„ j3 



Inserting this expression into eq.( p.3|) and using eq.( p.4|) yields 



OO -I 

z = y 

^ ATI 

iV=0 ^ ^ • 



3/2' 



2-^(5) 



N 



•' 1=1 



m \ 
2^) 



3/2 



N 



N 



(3.7) 
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It is convenient to introduce the dimensionless field 



x) = 2mJnG/3 ^{x) . 



(3.8) 



Then, 



where 



Z 



V(j)e ^-ff 



eff 



An 



T 



(3.9) 



(3.10) 



The partition function for the gas of particles in gravitational interaction has been trans- 
formed into the partition function for a single scalar field 4>{x) with local action 



sm-)] 



eff 



(TX 



(3.11) 



The (f) field exhibits an exponential self- interaction — /i^ e"^*^^-*. 

Notice that the effective temperature Tejj for the 0-field partition function turns out to 
be inversely proportional to T whereas the characteristic length behaves as ~ 7^-1/4 _ 
This is a duality-type mapping between the two models. 

It must be noticed that the term — e'^'-^^ makes the 0-field energy density unbounded 
from below. Actually, the initial Hamiltonian (|3.1|) is also unbounded from below. This 
unboundness physically originates in the attractive character of the gravitational force. In- 
cluding a short-distance cutoff [see sec. 2A, below] eliminates the zero distance singularity 
and hence the possibility of zero-distance collapse which is unphysical in the present context. 
We therefore expect meaningful physical results in the cutoff theory. Moreover, assuming 
zero boundary conditions for {/)(r) at r ^ 00 shows that the derivatives of cj) must also be 
large if e*^ is large. Hence, the term i(V0)^ may stabilize the energy. 

The action (|3.11|) defines a non-renormalizable field theory for any number of dimensions 
D > 2. This is a further reason to keep the short- distance cutoff non-zero. 

Let us compute now the statistical average value of the density p(r) which in the grand 
canonical ensemble is given by 



00 1 

< p(f) >= z~' y: 



N=0 



m \ 
2^) 



\ 3/2' 



N 



N 



Yl d^qi pif) e 



<ljtj<N \qi-qj\ 



1 = 1 



(3.12) 



As usual in the functional integral calculations, it is convenient to introduce sources in 
the partition function (|3.9| ) in order to compute average values of fields 



— I d.^^[h(^^)^ ]+ / d^x J{x) 



(3.13) 



The average value of (f){r) then writes as 
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< r >= 



ShgZ 
6J{f) 



In order to compute < p(r) > it is useful to introduce 



V[J{.)]^l(3Gm' [ 



Then, we have 



p(r) e 



\x-y\>a \X - y\ 
1 _ 



[p{x) + J{x) ] [p{y) + J{y) 



Tl • X 

eff 



5J{f 



By following the same steps as in eqs.(^]6| 

< p(f) > -- 



we find 



m 



3/2 



2n(3 



N 



m 



-1 



e 



i (V5)2-2m^7rG/3{(x) J{x) 



N 

f ... f \{d\i e^-v^EL 
•' 1=1 



1 vv2 / ^ 
y2 / 



rri • r 



5J{r) 



logZ[J(.)] 



I j=o 



Performing the derivatives in the last formula yields 



< pif) >= -— — l^'f^ V^(j){r) e ^-s! •' l^y ' j 2;[o] 

Teff J J 

We then see that in the 0-field language the particle density expresses as 

1 



p(f) 



eff 



By using eq.( p.l8|) , the gravitational potential at the point r 

d^x 



U{r) = -Gm J 



X — r\ 



can be expressed as 



U(r) = -— 4>(r) . 
m 



We can analogously express the correlation functions as 
C(ri, f2) = < p(ri)p(f2) > - < p(fi) >< p(f2) > 



\Teff 

This can be also written as 



I j=o 
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C{ri,r2) 



fx" 



-eff 



(3.21) 



A simple short distance regularization of the Newtonian force for the two-body potential 



IS 



Va(r 



[l-e{a-r)] 



6{x) being the step function. The cutoff a can be chosen of the order of atomic distances 
but its actual value is unessential. 

The introduction of the short-distance cutoff elliminates the unphysical short distance 
collapse of the gravitational gas. please notice, that the interaction between atoms or 
molecules is repulsive and not gravitational for short distances (Van der Waals forces, for 
instance). 



A. D-dimensional generalization 



This approach generalizes immediately to D-dimensional space where the Hamiltonian 
12]) takes then the form 

N 2 1 

HN = j:^-Gm' Y: T^^m^^ iovD^2 (3.22) 



and 



^ 7)2 1 

HN = y —-Gm^ y log— — , atD = 2. (3.23) 



The steps from eq.( |3.1| ) to ( p.9| ) can be trivially generalized with the help of the relation 

1 



V V \d 2 = 5{x - y) (3.24) 
F — y| 

in Z)-dimensions and 

log — — - = -C2 5{x - y) 

\x - y\ 

aX D = 2. 
Here, 

Cd = {D~ 2) — ^ for Dy^2 and C2 = 2n . (3.25) 



We finally obtain as a generalization of eq.( |3.9| ). 



Z^l I P^e-'*?/'"'-!*'™'-'''"""!, (3.26) 



where 
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We have then transformed the partition function for the L>-dimensional gas of particles 
in gravitational interaction into the partition function for a scalar field (p with exponential 
interaction. The effective temperature T^jf for the 0-field partition function is inversely 
proportional to T for any space dimension. The characteristic length behaves as ~ 

rp-{D-2)/4_ 



IV. SCALING BEHAVIOUR 

We derive here the scaling behaviour of the field following the general renormalization 



group arguments in the theory of critical phenomena |T5|,|T6| 



A. Classical Scale Invar iance 



Let us investigate how the action transforms under scale transformations 

X ^ xx = \x , (4.1) 



where A is an arbitrary real number. 

In D-dimensions the action takes the form 



Sim /rf^x[i(V0)^ -/.2e^(^^) 

^eff J 



(4.2) 



.2 

We define the scale transformed field (t)\{x) as follows 

(t)x{x)=(t){\x) + \ogW (4.3) 

Hence, 

(V0a(^))' = (V.,0(fA))' , e^^^^^ = \^ e^^'^^ 
We find upon changing the integration variable in eq.( [4.2|) from x to x\ 

S[U-)] = A^-^ ^[0(0] (4.4) 

We thus see that the action ( [4. 21) scales under dilatations in spite of the fact that it 
contains the dimensionful parameter /i^. This remarkable scaling property is of course a 
consequence of the scale behaviour of the gravitational interaction in D dimensions 0. 

In particular, in D = 2 the action ( [4. 21) is scale invariant. In such special case, it is 
moreover conformal invariant. 



The (Noether) current associated to the scale transformations (|4.1| ) is 

Ji{x) = X, Tij{x) + 2 Vi(j){x) , (4.5) 
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where Tij{x) is the stress tensor 



Tij{x) = Vi(f){x) Vj(f){x) - 6ij L 



and L = |(V0)^ — /i^ e"^^^^ stands for the action density. That i 



IS, 



Mx) = (x.V0 + 2) VMx) - X, [2 (V0)' - fi' e^(^^) 

By using the classical equation of motion we then find 

ViJiix) = (2 - D)L . 

This non-zero divergence is due to the variation of the action under dilatations [eq. ([4.4| )]. 
If 0(x) is a stationary point of the action ( [4 .21 ): 

VV(^) + /i' e"^^"^ = , (4.6) 
then (p\{x) [defined by eq.( |4.3|) ] is also a stationary point: 

VVa(x) e'^^(^) = . 

A rotationally invariant stationary point is given by 

2{D-2) 



(r) = log ■ 



(4.7) 



This singular solution is invariant under the scale transformations ([4.3|) . That is 

</>a(0 = ^\r) . 

Eq.( [4.7| ) is dilatation and rotation invariant. It provides the most symmetric stationary 
point of the action. Notice that there are no constant stationary solutions besides the 
singular solution 0o = —00. 



B. Thermal Fluctuations 



In this section we compute the partition function eqs. (|3.9| ) and (|3.13|) by saddle point 
methods. 

Eq. (|4.6| ) admits only one constant stationary solution 



bQ = -00 



(4.8) 



In order to make such solution finite we now introduce a regularization term e fi"^ 0(x) 
with e << 1 in the action S [eq.( |3.1lD ]. This corresponds to an action density 



L = 2(V0)' + w(0) 



(4.9) 
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where 



This extra term can be obtained by adding a small constant term — e fi'^/Tf.ff to p{x) in 
eqs.(^.4[) - (^.61). This a simple way to make 0o finite. 

We get in this way a constant stationary point at 0o = loge where m'(0o) = 0. However, 
scale invariance is broken since u"{(j)o) = — e /i^ 7^ 0. We can add a second regularization 



term to ^ 5 /i^ 



[X 



to L, (with 6 << 1) in order to enforce m"(0o) = 0. This quadratic 



term amounts to a long-range shielding of the gravitational force. We finally set: 

, 1 



[X 



6 (f){xy 



U{(f)) = -/i- 

where the two regularization parameters e and 6 are related by 

e{6)=6[l-\ogS], 
and the stationary point has the value 

00 = log (5 . 

Expanding around 0o 

0(f) = 00 + 6- xix) 



where g = J '^'^eff and x(f) is the fiuctuation field, yields 



1 



1 



(Vx)' 



9' 2 ■ • 
We see perturbatively in g that x(f) is a massless field. 



1 
2 



(4.10) 



Concerning the boundary conditions, for the ISM we must consider the system inside 
a large sphere of radius R (10~^ — 10~^ pc < R < 100 pc). That is, all integrals are 
computed over such large sphere. 



Using eq. ( p.l8| ) the particle density takes now the form 



p(f) 



1 



9 



o9Xif) 



9Xir} 



Teff Teff Teff 

It is convenient to renormalize the particle density by its stationary value 6 = e'^ 



p{r) 



9' 



D 



9X{r) 



(4.11) 



We see that in the S ^ limit the interaction in eq.( 4.1U| ) vanishes. No infrared divergences 
appear in the Feynman graphs calculations, since we work on a very large but finite volume 
of size R. Hence, in the S —>■ limit, the whole perturbation series around 0o reduces to the 
zeroth order term. 
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The constant saddle point 0o fails to catch up the whole field theory content. In fact, 
more information arises perturbing around the stationary point 0'^(r) given by eq.{L7_) |23 . 

Using eqs.( |3.21| ), ( ^.10| ) and ( [4.11| ) we obtain for the density correlator in the 6^0 
limit, 



C(ri,f2) 



2D 



9" 



exp 



Cd (/i|ri-f2|) 



D-2 



Cd (/i In - r2|) 



D-2 



For large distances, we find 



2 CI In - f2|2{^-2) 



H-0( 



ri - r2 



-3(D-2) 



(4.12) 



That is, the 0-field theory scales. Namely, the theory behaves critically for a continuum 
set of values of n and Tg//. 

Notice that the density correlator C{fi,f2) behaves for large distances as the correlator 
of This stems from the fact that is the most relevant operator in the series 

expansion of the density ( |4.11D 



(4.13) 



As remarked above, the constant stationary point 0o = log*^ ~^ ^oo only produces 
the zeroth order of perturbation theory. More information arises perturbing around the 
stationary point 0'^(r) given by eq.( [4.7| ) [ p3[ . 

We succeded in this way to to construct a scale invariant correlations in perturbation 
theory for the expenential interaction. The existence of such expansion is a necessary con- 
dition to show that the theory is scale invariant. Perturbation theory around the stationary 
point 



) shares such scale invariant property |23 



C. Renormalization Group and Finite Size Scaling Analysis 

As is well known [pTSH TTH , physical magnitudes for infinite volume systems diverge at the 
critical point as A to a negative power. A measures the distance to the critical point. (In 
condensed matter and spin systems, A is proportional to the temperature minus the critical 
temperature [p!6ip!7|). One has for the correlation length ^, 



and for the specific heat C, 

C ~ A-" . (4.14) 

Correlation functions scale at criticality. For example, the scalar field (which in spin 
systems describes the magnetization). 



< 0(r)0(O) >~ r 



-l—rj 
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The critical exponents z/, a and t] are pure numbers that depend only on the universality 
class [llMl. 



For a finite volume system, all physical magnitudes are finite at the critical point. 
Indeed, for a system whose size R is large, the physical magnitudes take large values at the 
critical point. Thus, for large R, one can use the infinite volume theory to treat finite size 
systems at criticality. In particular, the correlation length provides the relevant physical 
length ^ ~ i?. This implies that 

A ~ R-^/" . (4.15) 

We can apply these concepts to the 0-theory since, as we have seen in the previous section, 
it exhibits scaling in a finite volume ~ R^. Namely, the two points correlation function 
exhibits a power-like behaviour in perturbation theory as shown by eq. ([4.12|) . This happens 
for a continuum set of values of T^jf and fi^. Therefore, changing ^'^/T^ff keeps the theory 
in the scaling region. At the point ^"^/T^ff = 0, the partition function Z is singular. From 
eq. (|3.10|) , we shall thus identify 



Teff V 277 / 

Notice that the critical point A = 0, corresponds to zero fugacity. [For simplicity, we restrict 
from now on to the three dimensional case]. 

Thus, the partition function in the scaling regime can be written as 

Z(A)= / / P0e-^*+^/'^'^^'''' , (4.17) 



where S* stands for the action ( p.ll| ) at the critical point A = 0. 
We define the renormalized mass density as 

m p{x)ren = m e^^"^^ (4.18) 

and we identify it with the energy density in the renormalization group. [Also called the 
'thermal perturbation operator'.] This identification follows from the fact that they are the 
most relevant positive definite operators. Moreover, such identification is supported by the 
perturbative result ( |4.13| ). 



In the scaling regime we have [16| for the logarithm of the partition function 



^ log^(A) = j-^ r A'-" + F{A) , (4.19) 

V (2 — a;)(l — a) 

where F{A) is an analytic function of A around the origin 

F(A) =Fo + aA+^bA^ + ... . 

V = R^ stands for the volume and Fq, K, a and b are constants. 

Calculating the logarithmic derivative of Z{A) with respect to A from eqs.( [4.17D and 
from (|4.19|) and equating the results yields 
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(4.20) 



where we used the scahng relation a = 2 — vD [0,0] • 

We can apply here finite size scaling arguments and replace A by ~ R~~ [eg . (14.151 )], 

OA 1 — a 

Recalling eq.( |4.18| ), we can express the mass contained in a region of size R as 

i-R 



M{R) =mj e'^(^) d^x . 



Using eq.( |4.20|) we find 



< M{R) >=mVa + m 



K 



1 — a 



R^ + ... . 



and 



< p(r) >= ma + m 



K 



z/(l — a)f2o 



f V 



for r of order ~ R. 



(4.21) 



(4.22) 



where fi^ is the surface of the unit sphere in D-dimensions. 

The energy density correlators obey the renormalization group equations at criticality 
||T6| - Ijl^. One has for the connected L-point correlations of the energy density ('thermal 
operator'). 



^ d 1 



.k=l 



< p{xi) . ..p{xl) >' 







(4.23) 



for the three dimensional space. The homogeneity assumption implies that < p{xi)p{x2) > 
is only a function of the difference Xi — X2 for large jxi — X2|. We can therefore write for the 
density-density correlators (|3.2CI| ) in 3 space dimensions using by eq. ([4.23|) for L = 2, 



C(fi,f2) ~ |fi - ral- 



^-6 



(4.24) 



where both fi and r2 are inside the finite volume ~ R^. 

The perturbative calculation (|4.12|) matches with this result for = |- That is, the 
mean field value for the exponent u. 

Let us now compute the second derivative of log 2(A) with respect to A in two ways. 
We find from eq. ( |4.19| ) 



dA- 



logZ(A) = V A-'^K + b 



We get from eq. ( |4.17| ) 



9A2 



log Z{A) = J d^x d^y C{x, y)^ R^ J 



R d^ 



X 



X 



2D-2dH 



A-^ ~ R^ A-" 



(4.25) 
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where we used eq. ( [4. 151) , eq. (f4.24| ) and the scahng relation a = 2 — vD [lT6| , p^ . We conclude 
that the scaling behaviours, eq.( 4.19| ) for the partition function, eq.( 4.14 ) for the specific 
heat and eq. ( [4.24| ) for the two points correlator are consistent. In addition, eqs.( [4.21| ) and 
( [4.25| ) yield for the mass fluctuations squared 



(AM(i?))2 = < > - < Af >2 ~ j d^x d'^y C{x, y) ~ R^'^" 



Hence, 



AM(i?) ~ R^" 



(4.26) 



The scaling exponent v can be identified with the inverse Haussdorf (fractal) dimension 
dfj of the system 

dH = - ■ 

V 

In this way, M ~ R^^ according to the usual definition of fractal dimensions |T^ . 
The Fourier transform, 

p(fc) = y d^r e^-'' p{f) , 

will scale as: 

< p(fc) >~ k--^ . (4.27) 
This quantity is related to the number of clouds with mass M , n{M) , as follows 



n(M) ~ < p(M) > ~ M" 
From eqs.( p.lOI ) and ( [4.28| ) we find 

X = dn , 

which is in agreement with H. 



(4.28) 



Using eq.( [4.24| ) we can compute the average potential energy as 

d^x d^y 



< V >=-pGm^ f 



Cix,y)r^R--' 



i\x-y\>a \X — y\ 

From here and eq. (|4.26| ) we get as virial estimate for the atoms kinetic energy 



2 <^> oi-i 
< V >= ~ R" 

M{R) ■ 
This corresponds to a velocity dispersion 

That is, we predict [see eq. (|2.1|)] a new scaling relation 

1/1 \ 1. , 



(4.29) 



(4.30) 



The calculation of the critical amplitudes [that is, the coefficients in front of the powers 
of R in eqs. (|4.26| ), ( [4.24|) and (|4.3CI|) ] is beyond the scope of the present paper |l23[ . 



The 0-field connected correlators obey at criticality the renormalization group equation 
111-0 



N 



N 



Lfe=l 



< 0(fi)...0(fi) >' 



0, 



(4.31) 



for the three dimensional space. 



It follows from eq. ([4.31|) that the two-point subtracted correlator of the gravitational 
potential behaves for long distances as 



< (f){xi)(f){x2) > ~ |xi - X2I 



-1-rj 



D. Values of the scaling exponents and the fractal dimensions 

The scaling exponents u, a considered in sec IIIC can be computed through the renor- 
malization group approach. The case of a single component scalar field has been extensively 
studied in the literature |TB| , p7| , pil| . Very probably, there is an unique, infrared stable fixed 
point in three space dimensions: the Ising model fixed point. Such non-perturbative fixed 
point is reached in the long scale regime independently of the initial shape of the interaction 
[eq.Q] g. 

The numerical values of the scaling exponents associated to the Ising model fixed point 



are 



ly = 0.631... 



dn = 1.585... and a = 0.107... . 



(4.32) 



In the (j) field model there are two dimensionful parameters: fi and Tg//. The dimension- 
less combination 

g = fiTeff = (Stt) ' VZ j;^^ 

acts as the coupling constant for the non-linear fluctuations of the field (j). 

Let us consider a gas formed by neutral hydrogen at thermal equilibrium with the cosmic 
microwave background. We set T = 2.73 K and estimate the fugacity z using the ideal gas 
value 

^mT ) ^ 

Here we use p = 5q atoms cm~^ for the ISM density and 5q ~ 10^°. Eq. eqs. (|3.10|) yields 



- = 2.7 



AU~30AU and = pT, 



eff 



4.9 10 



-58 



510 



-53 



(4.33) 



This extremely low value for g^ suggests that the perturbative calculation [sec. IIIB] 
may apply here yielding the mean field values for the exponents, i. e. 



u = 1/2 



d 



H 



rj = and a = . 



(4.34) 
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That is, the effective couphng constant grows with the scale according to the renormahzation 
group ffow (towards the Ising fixed point). Now, if the extremely low value of the initial 
coupling eq. ( |4.33| ) applies, the perturbative result (mean field) will hold for many scales (the 
effective g grows roughly as the length). 

fi~^ indicates the order of the smallest distance where the scaling regime applies. A safe 
lower bound supported by observations is around 20 AU ~ 3. 10^^ cm , in agreement with 
our estimate. 



Our theoretical predictions for AM{R) and Av [eqs.( [4.26|) and (|4.30|) ] both for the Ising 
eq.( [4.32 ) and for the mean field values eq.( |4.34 ), are in agreement with the astronomical 
observations [eq. ( |2.1| )]. The present observational bounds on the data are larger than the 
difference between the mean field and Ising values of the exponents dn and q. 

Further theoretical work in the 0-theory will determine whether the scaling behaviour is 



given by the mean field or by the Ising fixed point 



E. The two dimensional gas and random surfaces fractal dimensions 



In the two dimensional case [D = 2) the partition function (|3.26|) describes the Liouville 
model that arises in string theory |^ and in the theory of random surfaces (also called two- 
dimensional quantum gravity). For strings in c-dimensional Euclidean space the partition 
function takes the form lESl 



(4.35) 



This coincides with eq.( p.26| ) aX D = 2 provided we flip the sign of /x^ and identify the 
parameters (|3.27|) as follows. 



12 ' ^ 



(4.36) 



Ref. states that = A for c < 1, = 3 for c = 2 and djj = 2 for c > 4. In our 
context this means 

2^ 11 

dH = 2 ioT T< — Gm^ , = 3 for T = 2 Gm^ and = 4 for T> — Gm^ . 

- 12 ' - 6 



For c — i> oo, ^ and we can use the perturbative result (|4.12| ) yielding = |, dn = 2 
in agreement with the above discussion for c > 4. 



F. Stationary points, the Poisson equation and the Jeans length 



The stationary points of the 0-field partition function ( ^.9] ) are given by the non-linear 
partial differential equation 



In terms of the gravitational potential U{x) [see eq. (|3.19| )], this takes the form 
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27T 



■U{r) 



(4.37) 



This corresponds to the Poisson equation for a thermal matter distribution fulfilling an ideal 
gas in hydrostatic equilibrium, as can be seen as follows |T^. The hydrostatic equilibrium 
condition 

VP(f) = -m p(f) Vf/(f) , 
where P(r) stands for the pressure, combined with the equation of state for the ideal gas 

P = Tp, 

yields for the particle density 

p(f) = Po e'"^^^^ , 
where po is a constant. Inserting this relation into the Poisson equation 

V^f/(f) = A-nGmp{r) 

yields eg. ( [4.37|) with 



Po 



For large r, eq.( |4.37D gives a density decaying as 



p(f) 



T 1 



2-KGm 



O 



U{r) 



T 



m 



log 



2nGpo 2 
— 7^ — r 
T 



+ 



(4.38) 



. (4.39) 



Notice that this density, which describes a single stationary solution, decays for large r 
faster than the density (f4.22|) governed by thermal fiuctuations. 

Spherically symmetric solutions of eq. (|4.37|) has been studied in detail The small 
fiuctuations around such isothermal spherical solutions as well as the stability problem were 



studied in 22 



The Jeans distance is in this context, 

1 _ vg(2.)^/' ^ 

m ^/GW^o ^/iGmy^T'^l^ 
This distance precisely coincides with [see eq.( p.l(]| )] up to an inessential numerical 
12/7r). Hence, p, the only dimensionful parameter in the 0-theory can be 




coefficient 

interpreted as the inverse of the Jeans distance. 

We want to notice that in the critical regime, dj grows as 

since po = A ~ R~'^" vanishes as can be seen from eqs.( |4.15| ), ( [4.16|) and ( |4.38| ). In this 
estimate we should use for consistency the mean field value dn = 2, which yields dj ~ R. 

This shows that the Jeans distance is of the order of the size of the system. The Jeans 
distance scales and the instability is therefore present for all sizes R. 

Had dj being of order larger than R, the Jeans instability would be absent. 

The fact that the Jeans instability is present precisely at dj R is probably essential 
to the scaling regime and to the self-similar (fractal) structure of the gravitational gas. 
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V. GALAXY DISTRIBUTIONS 



One obvious feature of galaxy and cluster distributions in the sky is their hierarchical 
property: galaxies gather in groups, that are embedded in clusters, then in superclusters, 
and so on. (Shapley 1934, Abell 1958). Moreover, galaxies and clusters appear to obey 
scaling properties, such as the power-law of the two point-correlation function: 



^(r) oc r~ 



-7 



with the slope 7, the same for galaxies and clusters, of ~ 1.7 (e.g. ^^). This scale-invariance 
has suggested very early the idea of fractal models for the clustering hierachy of galaxies (de 
Vaucouleurs 1960, 1970; Mandelbrot 1975). Since then, many authors have shown that a 
fractal distribution indeed reproduces quite well the aspect of galaxy catalogs, for example 
by simulating a fractal and observing it, as with a telescope (Scott, Shane & Swanson, 
1954; Soneira & Peebles 1978). Sometimes the analysis has been done in terms of a multi- 
fractal medium (Balian & Schaeffer 1989, Castagnoli & Provenzale 1991, Martinez et al 
1993, Dubrulle & Lachieze-Rey 1994). 

There is some ambiguity in the definition of the two-point correlation function ^(r) 
above, since it depends on the assumed scale beyond which the universe is homogeneous; 
indeed it includes a normalisation by the average density of the universe, which, if the 
homogeneity scale is not reached, depends on the size of the galaxy sample. Once ^(r) 
is defined, one can always determine a length tq where ^(/"o) =1 (Davis & Peebles 1983, 
Hamilton 1993). For galaxies, the most frequently reported value is tq ^ 5h~^ Mpc (where 
h = i^o/lOOkm s^^Mpc^^), but it has been shown to increase with the distance limits of 
galaxy catalogs (Davis et al 1988). tq is called 'correlation length' in the galaxy literature. 
[The notion of correlation length ^0 is usually different in physics, where ^0 characterizes the 
exponential decay of correlations (~ e~^^^°). For power decaying correlations, it is said that 
the correlation length is infinite]. 

The same problem occurs for the two-point correlation function of galaxy clusters; the 
corresponding ^(r) has the same power law as galaxies, their length has been reported 
to be about ro ~ 25h~^ Mpc, and their correlation amplitude is therefore about 15 times 
higher than that of galaxies (Postman, Geller & Huchra 1986, Postman, Huchra & Geller 
1992). The latter is difficult to understand, unless there is a considerable difference between 
galaxies belonging to clusters and field galaxies (or morphological segregation). The other 
obvious explanation is that the normalizing average density of the universe was then chosen 
lower. 

This statistical analysis of the galaxy catalogs has been criticized by Pietronero (1987), 
Einasto (1989) and Coleman & Pietronero (1992), who stress the unconfortable dependence 
of ^ (r) and of the length tq upon the finite size of the catalogs, and on the a priori assumed 
value of the large-scale homogeneity cut-off. A way to circumvent these problems is to deal 
instead with the average density as a function of size (cf §2). It has been shown that the 
galaxy distribution behaves as a pure self-similar fractal over scales up to ~ 100h~^ Mpc, the 
deepest scale to which the data are statistically robust (Sylos Labini et al 1996; Sylos Labini 
& Pietronero 1996). This is more consistent with the observation of contrasted large-scale 
structures, such as superclusters, large voids or great walls of galaxies of ~ 200/i~^ Mpc 
(de Lapparent et al 1986, Geller & Huchra 1989). After a proper statistical analysis of all 
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available catalogs (CfA, SSRS, IRAS, APM, LEDA, etc.. for galaxies, and Abell and AGO 
for clusters) Pietronero et al (1997) state that the transition to homogeneity might not yet 
have been reached up to the deepest scales probed until now. At best, this point is quite 
controversial, and the large-scale homogeneity transition is not yet well known. 

Isotropy and homogeneity are expected at very large scales from the Cosmological Prin- 
ciple (e.g. Peebles 1993). However, this does not imply local or mid-scale homogeneity 
(e.g. Mandelbrot 1982, Sylos Labini 1994): a fractal structure can be locally isotropic, but 
inhomogeneous. The main observational evidence in favor of the Cosmological Principle is 
the remarkable isotropy of the cosmic background radiation (e.g. Smoot et al 1992), that 
provides information about the Universe at the matter/radiation decoupling. There must 
therefore exist a transition between the small-scale fractality to large-scale homogeneity. 
This transition is certainly smooth, and might correspond to the transition from linear per- 
turbations to the non-linear gravitational collapse of structures. The present catalogs do 
not yet see the transition since they do not look up sufficiently back in time. It can be 
noticed that some recent surveys begin to see a different power-law behavior at large scales 
(A ^ 200 - 400/i"i Mpc, e.g. Lin et al 1996). 

There are several approaches to understand non-linear clustering, and therefore the dis- 
tribution of galaxies, in an infinite gravitating system. Numerical simulations have been 
widely used, in the hope to trace back from the observations the initial mass spectrum of 
fiuctuations, and to test postulated cosmologies such as CDM and related variants (cf Os- 
triker 1993). This approach has not yet yielded definite results, especially since the physics 
of the multiple-phase universe is not well known. Also numerical limitations (restricted 
dynamical range due to the softening and limited volume) have often masked the expected 
self-similar behavior (Colombi et al 1996). A second approach, which should work essen- 
tially in the linear (or weakly non-linear) regime, is to solve the BBGKY hierarchy through 
closure assumptions (Davis & Peebles 1977; Balian & Schaeffer 1989). The main assump- 
tion is that the A^-points correlation functions are scale-invariant and behave as power-laws 
like is observed for the few-body correlation functions. Crucial to this approach is the de- 
termination of the void probability, which is a series expansion of the A^-points correlation 
functions (White 1979). The hierachical solutions found in this frame agree well with the 
simulations, and with the fractal structure of the universe at small-scales (Balian & Schaef- 
fer 1988). A third approach is the thermodynamics of gravitating systems, developped by 
Saslaw & Hamilton (1984), which assumes quasi thermodynamic equilibrium. The latter is 
justified at the small-scales of non-linear clustering, since the expansion time-scale is slow 
with respect to local relaxation times. Indeed the main effect of expansion is to subtract the 
mean gravitational field, which is negligible for structures of mean densities several orders 
of magnitude above average. The predictions of the thermodynamical theory have been 
successfully compared with N-body simulations (Itoh et al 1993), but a special physical pa- 
rameter (the ratio of gravitational correlation energy to thermal energy) had to be adjusted 
for a better fit (Bouchet et al 1991, Sheth & Saslaw 1996, Saslaw & Fang 1996). 

We present in a new approach based on field theory and the renormalisation group 
to understand the clustering behaviour of a self-gravitating expanding universe. We also 
consider the thermodynamics properties of the system, assuming quasi-equilibrium for the 
range of scales concerned with the non-linear regime and virialisation. We find an exact 
mapping between the self-gravitating gas and a continuous field theory for a single scalar 
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field with an exponential self-coupling. This allows us to use statistical field theory and the 
renormalisation group to determine the scaling behaviour. The small-scale fractal universe 
can be considered critical with large density fluctuations developing at any scale. We derive 
the corresponding critical exponents. They are very close to those measured on galaxy 
catalogs through statistical methods based on the average density as function of size; these 
methods reveal in particular a fractal dimension D ^ 1.5 — 2 (Di Nella et al 1996, Sylos 
Labini & Amendola 1996, Sylos Labini et al 1996). This fractal dimension is strikingly 
close to that observed for the interstellar medium or ISM (e.g. Larson 1981, Falgarone 
et al 1991). We show in ref. that the theoretical framework based on self-gravity that 
we have developped for the ISM (de Vega, Sanchez & Combes 1996a,b, hereafter dVSC) 
is also a dynamical mechanism leading to the small scale fractal structure of the universe. 
This theory is powerfully predictive without any free parameter. It allows to compute the 
A^-points density correlations without any extra assumption. 



VI. GALAXY CORRELATION FUNCTIONS AND MASS DENSITY IN A 

FRACTAL 



The use of the two point correlation function ^(r) widely spread in galaxy distributions 
studies, is based on the assumption that the Universe reaches homogeneity on a scale smaller 
than the sample size. It has been shown by Coleman, Pietronero & Sanders (1988) and 
Coleman & Pietronero (1992) that such an hypothesis could perturb significantly the results. 
The correlation function is defined as 

_ < n{ri).n{ri + r) > _ ^ 
< n >2 

where n[r) is the number density of galaxies, and < ... > is the volume average (over d^Vi). 
The length ro is defined by ^(^o) = 1. The function ^(r) has a power- law behaviour of slope 
—7 for r < ro, then it turns down to zero rather quickly at the statitistical limit of the 
sample. This rapid fall leads to an over-estimate of the small-scale 7. Pietronero (1987) 
introduces the conditional density 

^ _ < n{ri).n{ri + r) > 
<n> 

which is the average density around an occupied point. For a fractal medium, where the 
mass depends on the size as 

M(r) oc 

D being the fractal (Haussdorf) dimension, the conditional density behaves as 

r(r) oc r^-^ 

This is exactly the statistical analysis used for the interstellar clouds, since the ISM as- 
tronomers have not adopted from the start any large-scale homogeneity assumption (cf 
Pfenniger & Combes 1994). 
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The fact that for a fractal the correlation ^(r) can be highly misleading is readily seen 
since 



and for a fractal structure the average density of the sample < n > is a decreasing function 
of the sample length scale. In the general use of ^(r), < n > is taken for a constant, and we 
can see that 

D = 3-7 . 

If for very small scales, both ^(r) and r(r) have the same power-law behaviour, with the 
same slope —7, then the slope appears to steepen for ^(r) when approaching the length tq. 
This explains why with a correct statistical analysis (Di Nella et al 1996, Sylos Labini & 
Amendola 1996, Sylos Labini et al 1996), the actual 7 ^ 1 — 1.5 is smaller than that obtained 
using C,{r). This also explains why the amplitude of ^(r) and tq increases with the sample 
size, and for clusters as well. 

In the following, we adopt the framework of the fractal medium that we used for the 
ISM (dVSC), and will not consider any longer ^(r). 



VII. EQUATIONS IN THE COMOVING FRAME 

Let us consider the universe in expansion with the characteristic scale factor a{t). For 
the sake of simplicity, we modelise the galaxies by points of equal masses m, although they 
have a mass spectrum (it may be responsible for a multi-fractal structure, see Sylos Labini 
& Pietronero 1996). 

The present analysis can be generalised to galaxies of different masses following the lines 
of sec. IX 0. We expect to come to this point in future work. 

If the physical coordinates of the particles are r, we can introduce the comoving coordi- 
nates X such that 

r = a{t) X 

The Lagrangian for a system of particles interacting only by their self-gravity can be 
written as 



N 



ln=y: 



i=l 



a{t) Xi-—- (j){xi{t)) 
2 a[t) 



(7.1) 



where 0(x) is the gravitational potential in the comoving frame, determined by the Poisson 
equation 

V20(x) =47rG'p(x,t) , (7.2) 
and p(x, t) is the mass density. For our system of point particles, 

N 

p{x,t) = m '^6{x — Xi(t)) (7.3) 

i=l 

and therefore the solution of the Poisson equation takes the form 
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N ^ 

(x) = -Gm y - — —-r . (7.4) 



The canonical momenta and Hamiltonian of the system are 

Pi = m a{tY Xi 



N 



hn=y: 



i=l 



1 _ 9 m 



Pi + —r?: 'Pi.^i^^)) 



2ma{tf a{t) 



2ma(t)2 a{t) ^^i^^^ ' 

We see that the A^-particle Hamiltonian in cosmological spacetime eq.( [7.5|) can be obtained 
from the Minkowski Hamiltonian [a{t) = 1] by making the replacements 

m-^ma{tf , G ^ G a{t)-\ (7.6) 

As a first approximation, we shall assume in the following that the characteristic time 
of the particle motions under the gravitational self-interaction are shorter than the time 
variation of a{t). We can then consider that this system of self-gravitating particles is at 
any time in approximate thermal equilibrium. This hypothesis is true of course for structures 
that have already decoupled from the expansion, and are truly self-gravitating and virialised. 
It could be also valid for the whole non-linear regime of the gravitational collapse. As for 
the linear regime, we know already that the primordial fluctuations are not forgotten in 
the large-scale structures, and therefore the resulting correlations will depend on initial 
conditions, and not be entirely determined by self-gravity. 

The above assumption introduces a natural upper limit in the scales concerned by the 
theory developped below. Similarly to the case of the interstellar medium, the fractal struc- 
ture considered is bounded by a short distance cut-off and by a large-scale limit as well 
(dVSC). 

The short distance cut-off corresponds to the appearence of other physics at short scale, 
essentially dissipative, which we do not need to introduce. In addition, the short distance 
cut-off avoids the gravo-thermal catastrophe. For the ISM, the cut-off was naturally the size 
of the smaller fragments, of the order of the Jeans length. Here the cut-off corresponds also 
to the size of the 'particles' considered, i.e. the galaxy size, below which another physics 
steps in, related to stellar formation and radiation. 

The fact that in the catalogs, we are observing in projection large-scale structures at 
different epochs, with different values of the scale factor a{t), could slightly modify the 
fractal dimension. Even though fractal structures are self-similar, and scale-independent, 
the largest scales are systematically observed at a younger epoch where the contrast has not 
grown up as high as today. This evolution effect however should be significant only at high 
redshift (> 1), and the present catalogs are not yet statistically robust so far back in time 
(the average redshift of optical catalogs is about 0.1). 
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VIII. APPLICATION OF RENORMALIZATION GROUP THEORY 



As in all scale-independent problems, where the fluctuations cannot be represented by 
analytical functions, the renormalization group theory developped in the 1970's for the study 
of critical phenomena, appears here perfectly adapted (e.g. ||l5l). We can consider the fractal 
structure of the Universe as the critical state of a system, where fluctuations develop at any 
scale, with a very large correlation length (asymptotically infinite). The fluctuations that 
are distributed as a fractal of dimension D are the large-scale structures of the universe (cf. 

ill). 

We generalise now the study of an A^-body system only interacting through their own 
self-gravity to the case where the bodies are on a cosmological background. 

Another approach has been proposed for galaxy correlations but it yields different 
critical exponents. 

Let us apply the theory to the system of galaxy points, already defined in the previous 
section. Since they are considered in approximate thermal equilibrium, we will use the grand 
canonical ensemble, that also allows a variable number of particles. The grand partition 
function of the system can be written as in eq. |3.1| ) where now is given by eq.( |7.5| ). 



The functional representation for the grand partition function can be easily generalized 
for an arbitrary scale factor a{t). After the changes specified above in eq.( |7.(j| ), the local 
action becomes 



i(V0)2-A(t)2 e^(^^) 



(8.1) 



Notice that all quantities depend on time through the scale factor a{t) only. There is no 
integration over t. 

The mass parameter n in the </)-theory gets effectively multiplied by the scale factor 
a{t). Since the Jeans length dj ^ fi~^ according to eg. (p. 14 ), in comoving coordinates dj 
effectively becomes 

^ 1 



d 



J 



n jj,a{t) 



as one could have expected. 

On the other hand, the dimensionless coupling constant 



is unchanged by the replacements of eq. ( |7.6| ). 

Therefore, for any fixed time t we find the same scaling behaviour, after making the 
replacement 

jji ^ H a{t) 

and keeping the coupling g unchanged. 

Thus, the renormalisation group and finite size scaling analysis of sees. IVC and IVD 
apply without essential changes to the galaxy distributions. Namely, the mass fluctuations 
AM(i?) inside a volume R, 



{AM{R)y = < > - < M >2 
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will scale as 



AM(i?) - i?^ 



(8.2) 



The scaling exponent v can then be identified as before with the inverse Haussdorf 
(fractal) dimension D of the system 




As usual in the theory of critical phenomena, there are only two independent critical 
exponents. All exponents can be expressed in terms of two of them: for instance the fractal 
dimension D = and the independent exponent rj, which usually governs the spin- 

spin correlation functions. The exponent rj appears here in the 0-field correlator (dVSC), 
describing the gravitational potential, that scales as 

< 0(r) > - r-^(i+'') 

The values of the critical exponents depend on the fixed point that governs the long range 
behaviour of the system. 

The value of the dimensionless coupling constant = fiTeff should decide whether the 
fixed point chosen by the system is the mean field (weak coupling) or the Ising one (strong 
coupling). At the tree level, we estimate g ~ where N is the number of points in a 
Jeans volume Sj. The coupling constant appears then of the order of 1, and we cannot 
settle this question without effective computations of the rcnormalisation group equations. 
At this point, the predicted fractal dimension D should be between 1.585 and 2. 



A. Three point and higher correlations 

Our approach allows to compute higher order correlators without any extra assumption. 
The two and three point densities, 

D(fi,r2) = < n(fi)n(f2) > 

-D(ri, f2, fa) = < n(fi) n(r2) n{fz) > , (8.3) 
can be expressed as follows in terms of the correlation functions: 
-D(fi,f2) =ni n2 + Ci2 
-D(fi, f2, fs) = ni 712 ^3 + rii C23 + na C12 + n2 C13 + C123 ■ (8.4) 

Here, 

Ui =< n{fi) > , i = 1,2,3 , 
and Cij and Cijk are the two and three point correlation functions, respectively, 

Cij = C{fi, fj) 
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The behaviour of n^, Cij and Cijk in the scahng regime follow from the renormalisation 
group equations at criticahty (de Vega, Sanchez & Combes, in preparation). If we do not 
impose homogeneity at all scales, we find. 



D-3 



< n{r) >~ A r 

^.Tx»r-2 2{Z)-3) 
6(ri, 72) C:^ B \ 

<^(n,r2,r3) ~ Cr/ ' , (8.5) 

where A, B and C are constants and D = \ jv. 

We can now derive the three point density behaviour when one point, say ri, is far away 
from the other two. We find from eqs.(|8.4j) and (^.5|), 



L)(ri,r2) ~ Ar^ n2 +Bri , 

^('"l,'"2,?^3) ^ (?^2 ?^3 + C'23) 

+ B r^(^'') (712 + m) + C rl^""-^^ (8.6) 

Notice that this expression is dominated by the first term since D — 3 < 0. 

Higher point distributions can be treated analogously in our approach. We find that the 
dominant behaviour in the A^-points density is 

^ ri»ri,2<i<N N{D-3) /o ^7^ 

C;(ri,r2, . . . ,rAr) ~ r^' ' (8.7) 

Notice that when homogeneity is assumed to hold over all scales, the critical behaviour of 
the iV-point correlation function involves a factor ^['^^^^^^^^^ ^ |49 . 



Eqs.( |8.6[ - ^.7D are qualitatively similar, although not identical, to the behaviour inferred 
assuming the factorized hierarchical Ansatz (fhA), (Balian & Schaeffer 1989). That is, 

D{n, f^y^^ = n'{l + b r^f') (l + rf^^) , (8.8) 

I^(ri,r2,r3)^''^ = # {l + 6 (rfa^^^ + rg'^ + 



,0-3 

23 



+ Q3 

ri>>r2,r'3 3 



n 



^D-3 I ^D-3 ^D-3 , ^D-3 ^D-3 
'12 '13 "1" ' 12 '23 "T ' 13 '23 

1 + h rg-' + 2 rf -3 (6 + rg-') + Q3 rf""-'^ 



where ri2 = \fi — and so on. h and Q3 are constants. Notice that in the factorized 
hierarchical Ansatz, the fractal dimension D is not predicted but it is a free parameter. 

We see that the dominant behaviours in eqs.( |8.6| ) and ( ^.8] ) are similar in case the scaling 
exponents D — 3 are the same. 
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IX. DISCUSSION 



In previous sections we ignored gravitational forces external to the gas like stars etc. 
Adding a fixed external mass density Pext{'>^ amounts to introduce an external source 

J(f) = -T^ff pextir) , 

in eq.( |3.13| ). Such term will obviously affect correlation functions, the mass density, etc. 
except when we look at the scaling behaviour which is governed by the critical point. That 
is, the values we find for the scaling exponents dn and q are stable under external pertur- 
bations. 

We considered all atoms with the same mass in the gravitational gas. It is easy to 
generalize the transformation into the 0-field presented in section II for a mixture of several 
kinds of atoms. Let us consider n species of atoms with masses m^, 1 < a < n. Repeating 
the steps from eg. to yields again a field theory with a single scalar field but the 

action now takes the form 



^eff 



(9.1) 



where 



and m is just a reference mass. 

Correlation functions, mass densities and other observables will obviously depend on the 
number of species, their masses and fugacities but it is easy to see that the fixed points and 
scaling exponents are exactly the same as for the 0-field theory ( |3.9| )-( prT0| ). 

We want to notice that there is an important difference between the behaviour of the 
gravitational gas and the spin models (and all other statistical models in the same uni- 
versality class). For the gravitational gas we find scaling behaviour for a full range of 
temperatures and couplings. For spin models scaling only appears at the critical value of 
the temperature. At the critical temperature the correlation length ^ is infinite and the 
theory is massless. For temperatures near the critical one, i. e. in the critical domain, ^ is 
finite (although very large compared with the lattice spacing) and the correlation functions 
decrease as ~ e~^^^ for large distances r. Fluctuations of the relevant operators support 
perturbations which can be interpreted as massive excitations. Such (massive) behaviour 
does not appear for the gravitational gas. The ISM correlators scale exhibiting power-law 
behaviour. This feature is connected with the scale invariant character of the Newtonian 
force and its infinite range. 

The hypothesis of strict thermal equilibrium does not apply to the ISM as a whole where 
temperatures range from 5 to 50 K and even 1000 K. However, since the scaling behaviour is 
independent of the temperature, it applies to each region of the ISM in thermal equilibrium. 
Therefore, our theory applies provided thermal equilibrium holds in regions or clouds. 
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We have developped here the theory of a gravitationally interacting ensemble of bodies 
at a fixed temperature. In a real situation like the ISM, gravitational perturbations from 
external masses, as well as other perturbations are present. We have shown that the scaling 
solution is stable with respect to the gravitational perturbations. It is well known that 
solutions based on a fixed point are generally quite robust. 

Our theory especially applies to the interstellar medium far from star forming regions, 
which can be locally far from thermal equilibrium, and where ionised gas at lO^K together 
with coronal gas at lO^K can coexist with the cold interstellar medium. In the outer parts 
of galaxies, devoid of star formation, the ideal isothermal conditions are met 0]. Inside the 
Galaxy, large regions satisfy also the near isothermal criterium, and these are precisely the 
regions where scaling laws are the best verified. Globally over the Galaxy, the fraction of the 
gas in the hot ionised phase represents a negligible mass, a few percents, although occupying 
a significant volume. Hence, this hot ionised gas is a perturbation which may not change 
the fixed point behaviour of the thermal self-gravitating gas. 



In ref. 1^31 ^ connection between a gravitational gas of galaxies in an expanding universe 



and the Ising model is conjectured. However, the unproven identification made in ref. P? 



of the mass density contrast with the Ising spin leads to scaling exponents different from 
ours. 



Concerning the galaxy distributions, different scaling behaviours show up depending 
whether the system is homogeneous or not at large distances. The homogeneity property 
brings extra information which is not contained in the fundamental gravitational interaction. 
In condensed matter systems, such homogeneity hypothesis is easily verified by experiments. 
The homogeneity of the universe at large scales is a much more controversial issue 

Under the homogeneity hypothesis we find for the galaxy-galaxy correlator defined as in 
1^, 



) ^ < p(n))p(ro + f)> 
< p >^ 

where 7 = 6 — 2dH- This gives 'ymean field = 2 and •y ising = 2.830 .... Such numbers should 
be compared with the customary value 7 = 1.8 obtained from the observations [^ . 



In the inhomogeneous regime we find for the galaxy-galaxy correlations considered in ref. 

D{r) =< p{fo)p{fo + r ) > ~ r ^ 

where F = 3 — dn- This gives T mean field = 1 and Tjsing = 1.415 .... Such numbers should 
be compared with the observed value for F ~ 1 [^] (obtained for r > 30 Mpc). 

In conclusion, our renormalization group results seem to fit better with the analysis of 
the observations made in ref. [0] than with the standard lore . However, our predictions 



do not disagree in a dramatic way with the standard analysis of the observations. In both 
cases, mean field exponents give the better agreement with observations. 
Nuevas conclusiones: 

The statistical analysis of the most recent galaxy catalogs, without the assumption of 
homogeneity at a scale smaller than the catalog depth, has determined that the universe 
has a fractal structure at least up to ~ 100/i~^ Mpc (Sylos Labini et al 1996). The analysis 
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in terms of conditional density has revealed that the fractal dimension is between D = 1.5 
and 2 (Di Nclla ct al 1996, Sylos Labini & Amendola 1996). We apply a theory that we 
have developped to explain the fractal structure of the interstellar medium (dVSC), which 
has the same dimension D. The physics is based on the self-gravitating interaction of an 
ensemble of particles, over scales hmited both at short and large distances. The short- 
distance cut-off is brought by other physical processes including dissipation. The long-range 
limit is fixed by the expansion time-scale. In-bctwccn, the system is assumed in approximate 
thermal equilibrium. The dynamical range of scales involved in this thermodynamic quasi- 
equilibrium is at present limited to 3-4 orders of magnitude, but will increase with time. 

The critical exponents found in the theory do not depend on the conditions at the cut- 
off, which determine only the amplitudes. The theory is based on the statistical study of 
the gravitational field: it is shown that the partition function of the N-body ensemble is 
equivalent to the partition function of a single scalar field, with a local action. This allows 
to use field theory methods and the renormalisation group to find the scaling behaviour. 
We find scaling behaviour for a full range of temperatures and couplings. The theory then 
predicts for the system a fractal dimension D — 1.585 for the Ising fixed point, or D = 2 in 
the case of the mean-field fixed point. Both are compatible with the available observations. 
The iV-points density correlators are predicted to scale with exponent {N — 1){D — 3) when 
ri >> Tj, 2 <i < N. That is, -{N - 1) for the mean field, or —1.415 {N - 1) for the Ising 
point. 

We predict in addition a critical exponent — 1(1 + ?]) for the gravitational potential: that 
is, —0.500 for mean field or —0.519 for the Ising fixed point. 
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